Magnon decay in gapped quantum spin systems 
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In the 0(3) <T- model description of gapped spin systems, 5 = 1 magnons can only decay into 
three lower energy magnons. We argue that the symmetry of the quantum spin Hamiltonian often 
allows decay into two magnons, and compute this decay rate in model systems. Two magnon decay 
is present in Haldane gap S = 1 spin chains, even though it cannot be induced by any allowed 
term written in powers and gradients of the cr-model field. We compare our results with recent 
measurements of Stone et al. on a two-dimensional spin system. 

PACS numbers: 75.10.Jm, 75.40.Gb 



Introduction.- Stone et ali> have recently studied 
the spectrum of the 5 = 1 magnon excitations of 
the gapped quasi-two-dimensional spin dimer compound 
(C4Hi2N2)Cu2Cl6 (piperazinium hexachlorodicuprate, 
PHCC). They observed that the magnons become un- 
stable when their energy crosses the lower boundary of 
the two-magnon continuum, and measured magnon life- 
times. These observations are interesting because, at first 
glance, they appear to be in conflict with the standard 
phenomenological description of low energy S = 1 exci- 
tations in a confining spin-gap systems. This model 2 3 
is expressed in terms of a field n(x, r), where a; is a 
d-dimensional spatial coordinate, r is imaginary time, 
and n is three-component vector in spin space. Using 
Landau-Ginzburg arguments, the effective action for the 
magnons, near the minimum of their dispersion, can be 
written in the form 

S n = j d d xdr{^-({d Tn f+c\d xn f) + ^( n 2 - a ) 2 } 

Harmonic quantum fluctuations of n about n = in a 
renormalized effective potential then constitute the triply 
degenerate magnon in the spin-gap state. In this picture, 
it is clear from the quartic non-linearity in S n that this 
magnon can only decay into three magnons, and decay 
into two magnons is prohibited. 

It is useful to recall the constraints imposed by sym- 
metry more carefully for the familiar case of a spin 5* 
antifcrromagnetic chain (d = 1). This model has spin ro- 
tation symmetry, time reversal symmetry (under which 
r — > — t and n — > — ra), reflection symmetry about sites 
(under which x — > —x and n — > n), and reflection sym- 
metry about centers of bonds (under which, in the con- 
tinuum limit, x — > —x and n — » — n). These symmetries 
are very restrictive and only allow a single term with 
three powers of n: the topological '0-term— 



S„ = 



i6 J dxdr n ■ (d x n x d T n) , 



(1) 



in the fixed length limit u — > 00 with a = 1 and 9 quan- 
tized at 9 = S/2. In this limit the integrand in Sg is 
quantized in integer multiples of 47r, and only then is 



terms as they act on the fully renormalized quasiparticles 
of the gapped antiferromagnet with integer S. As noted 
above, these can be considered to be harmonic fluctu- 
ations of a renormalized field n about n = 0, and the 
amplitude of this renormalized field is not constrained 
to be unity. For such a field n, it is easy to see that 
there is no term with three powers of n, and arbitrary 
spatial and temporal gradients which is invariant under 
the symmetries. 

Here we will examine the issue of two magnon decay us- 
ing a lattice formulation. Our analysis begins below with 
a microscopic Hamiltonian appropriate for PHCC. We 
will compute the lifetime due to decay into two magnons, 
and compare our results with the observations of Stone et 
alX The confrontation between theory and experiments 
places new constraints on the values of the microscopic 
exchange constants. We then show that similar processes 
also exist for a generic S = 1 spin chain, and are consis- 
tent with all symmetries of the lattice Hamiltonian; this 
is the case even though such processes do not appear in 
any perturbation of the continuum theory. 

Magnon lifetimes in PHCC- The Cu 2+ spins in PHCC 
form a lattice of dimers in the (ac) crystalline plane as 
shown in Fig. ^ We will denote the interdimer exchange 
interactions as Jg" , where a, a 1 — 1,2 label two spins 
inside a dimer and the vector 8 = (m, n) = ma + nc 
connects the dimer centers. For a description of such 
systems it is convenient to use the so-called bond-boson 
formalism 4 where spin-i operators Si, S2 at each dimer 
(see Fig. ^) are represented in terms of three hardcore 
bosons t a , a — (x, y, z) which correspond to three excited 
triplet states of the dimer: 

I, + s 1 



(2) 



If one assumes the isotropic Heisenberg coupling between 
dimers, the spin Hamiltonian takes the form H = Hq +2 
+H 3 + H 4 , where H — JoJ^r^ia^ra i s determined by 
the intradimer exchange Jq, and iz2,3.4 correspond to the 
interdimer interaction. (The index r here labels the sites 
of the lattice formed by dimers). 

The quadratic part of the interdimer Hamiltonian 



1 invariant under all symmetries noted above. H 2 = - ^ ^-T"'(4aVa + trat r > a +h.c), (3) 



However, we arc interested in the nature of the non-linear 



(rr') 



as well as the four-particle interaction H4, 

U - I V* T cS lit ft i + _ +t + t + + \ 
4 ~ 2 Jl rr'\ L ru L r'a l r{3 l r l (3 VaV|3VaVj3J i 

{rr') 

depend only on specific combinations r , of exchange 
couplings between the dimers located at sites r and r', 

J cS , = -(J n , + J 22 , — J 12 , - J 21 ,) (4) 

"r—r' 2^ r — r' 1 "r — r' <-> T _ r i u r _ T ij. y^j 

The three-magnon term H3 depends on different com- 
binations of exchange constants F 1 — r < that are nonzero 
only if the interdimer coupling is not symmetric with re- 
spect to the interchange of the spin indices 1 and 2: 

#3 = 5 £ ^7 i? r-r-'(4a4'/3*r-' 7 ~ ^'jl^r-y) +h.C. 
(rr') 

F r — r i ~(J r _ r i J r — r i J r — r i ~\~ J r — r >)- (5) 

This type of interaction arises generically in the theory; 
however, it is usually neglected in practical calculations. 

If one neglects magnon interaction completely, taking 
into account only the quadratic part of the Hamiltonian 
Hq + H 2l and uses this approximation to describe the 
experimentally observed magnon dispersion and to find 
out the microscopic exchange couplings, one can obtain 
estimates for the intradimer exchange constant Jo and 
for "effective couplings" J^_ r i only. The individual cou- 
plings Jr- r >, Jr- r ', etc. cannot be determined in this 
nonintcracting magnon approximation that is commonly 
used by experimentalists for the description of the dis- 
persion data (see, e.g., Refs.000). 

Since an analytical treatment of the full interacting 
Hamiltonian for realistic models does not seem feasi- 
ble, one has to resort to extensive numerical calculations 
to extract the coupling information from the dispersion 
data. Such calculations based on cluster expansions were 
actually performed for various spin dimer materials (SaLifl 
however, even in this approach the problem remains diffi- 
cult because for weak interdimer interactions the magnon 
dispersion is rather insensitive to the variations of the 
individual exchange couplings which keep the "effective" 
interactions J^ r , constant. 

The magnon dispersion in PHCC was studied exper- 
imentally in Ref. |5| the spin excitations were found to 
have a spectral gap A « 1 meV and a bandwidth of about 
1.8 meV. It was shown that the observed spectrum is well 
described within the noninteracting magnon approxima- 
tion, if one assumes the existence of effective interdimer 
interactions Jf with <5 = (1,0), (0,1), (1,1), (1,-1), 
(2, 0), and (0, 2), with J ~ 2.33 meV, J^ 0) ~ 0.44 meV, 

J ( c ff 1} ~ 0.23 meV, J^ 1} ~ J^_ 1} ~ -0.084 meV, 
J ( e 2 ff o) - -0.073 meV, and J ( c ff 2) ~ -0.047 meV. 11 In this 
approximation, the magnon dispersion takes the form 

e(fe) - { J 2 + 2J B fc } 1 / 2 , B k =J2jfcos(k-S) (6) 

s 
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FIG. 1: A schematic view of the coupled spin dimer structure 
in the crystallographic (ac) plane of PHCC, and the notation 
used. Two spins in each dimer are distinguished by the spin 
index a = 1, 2, and each dimer is labeled by the position of 
its center r. Exchange couplings between individual spins are 

then denoted as J^Zj. In PHCC, the most important ex- 
change links are given by 1 — r' = (1,0), (0, 1), (1, 1), (1, — 1), 
and possibly also (2,0) and (0,2) (in lattice units). 

and the quasiparticle operators b k are connected to the 
dimer triplet-creating operators t k through the stan- 
dard Bogoliubov transformation t k = u k b k + v k b_ k , 
where the coefficients u k — cosh a:*, , v k — sinhxfc, with 
tanh(2x fc ) = -B k /(J Q + B k ) and sign(z; fc ) = -sign(-Bfe). 

The four-magnon interaction term H4 leads in fact to 
a renormalization of the amplitudes Jg S in Hi and can 
be ignored as far as one uses the values of the effective 
couplings Jg S obtained from a comparison with the ex- 
perimental data, since the above renormalization is al- 
ready included in this way. The same argument applies 
to the hardcore constraint which, for low density of vir- 
tual triplet pairs, in the leading approximation can be 
shown just to renormalize the self-energy of a magnon. 12 

The three-magnon term H3 does not contribute to the 
dispersion in the first order of the perturbation theory, 
and its main physical effect is to produce an instability of 
magnons with energies above the two-particle continuum 
threshold. In terms of transformed quasiparticle opera- 
tors bk this interaction can be written as 

(xyz) 12 

M12 = ^ Fg|(ui +2 - Vl+2)(V 1U2 - V2U1) 
S 

x sin(fci • S + k 2 ■ S) (7) 

- (ui - wi)(ui +2 u2 + vi +2 v 2 ) sin(fci • S) 

+ (u 2 - « 2 )(ui+2Ui +ui +2 wi)sin(fe 2 ■$)}, 

where (xyz) denotes a cyclic summation over the triplet 
component indices, (1,2) = fci.2, and the remaining 
terms describing three-particle creation are omitted for 
clarity. The decay rate T(p) of a magnon with the wave 
vector p according to the Fermi golden rule is 

T(p) =2ttJ2 \M k ^ k \ 2 S(e(p) - e{k) - e(p - k)) . (8) 

k 

After performing the delta-function integration, the lat- 
ter expression is reduced to a one-dimensional quadrature 
which can be computed numerically. 
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bond-boson calculation 

• exp. data by Stone el al. 



0.1 0.2 0.3 

wave vector (h, -h) in (ac) plane [r.l.u.] 



FIG. 2: Magnon linewidth T(p) in PHCC as a function of 
the wave vector p — (h, —h) [r.l.u.] in the (ac) plane. The 
solid line displays the result of the bond-boson calculation 
(see text) and the symbols correspond to experimental data 
of Stone et ali 



In Ref. © the magnon dispersion in PHCC was ana- 
lyzed on the basis of the fourth-order cluster expansion 
calculation, and it was concluded that a good fit to the 
data of Ref. can be achieved even without including the 
links with 8 = (2, 0) or (0, 2) (for those latter links, the 
existence of exchange paths seems questionable). Since 
we are not taking into account higher-order effects, we 
will adopt the coupling values J$ s proposed in Ref. for 
the sake of consistency; one may argue that this is not 
going to introduce any serious problem since the role of 
J| ff in our calculation here is just to mimic the correct 
behavior of e(k) and the Bogoliubov factors Uk, Wfe- 

Another conclusion of Ref. |9j was that the links with 
S = (1,1) and (1,-1) are fully due to the "diagonal" 
couplings of the J 12 and J 21 type, respectively, while the 
other intcrdimcr links have symmetric couplings in the 
sense that J 12 = Jj 1 and Jg 1 = Jg 2 , so that all F$ = 
except ^(1,-1) ~ ~ 0.09 meV. 

We have fitted the experimental data- for T as a func- 
tion of the wave vector p in the (ac) plane, using the am- 
plitudes Fg as fitting parameters. If one assumes that all 
the individual exchange couplings are antiferromagnetic, 
a natural constraint arises that the absolute value of the 
amplitude Fg cannot exceed the absolute value of the 
corresponding Jg S . It turns out that the observed mag- 
nitude of the linewidth 1 (the maximum of T ~ 0.13 meV 
is reached at p — (0, 0)) cannot be achieved if one adopts 
the coupling pattern proposed in Ref. [{J the contribution 
of ^(1,1) and alone is definitely insufficient (with 

r(0) one order of magnitude below the observed value). 
The linewidth is rather insensitive to the value of F(p,i)'- 
even at -F(o,i) = J(oi) the corresponding contribution to 
T(p = 0) is of the order of 10~ 4 meV. Similarly, inclusion 
of longer links like F ( 2i0 ) and -F(o,2) does not help, since 
at their "full strength" (i.e., at F$ = Jg S ) they yield the 
contributions to T(0) of the order of several neV. The 
only remaining way to explain the observed values of T 



is to assume that -F(io) 7^ 0; indeed, as shown in Fig. [21 a 
reasonable fit of the experimental data is achieved with 

F(i,o) = - 7J ( i ff o)> F (i,i) = ~ J {t,i)> and = J (°i ff -i)- 

The shown fit is not unique: the signs of F$ were chosen 
to provide the maximum enhancement of the total am- 
plitude of the matrix element M12, and the amplitudes 
F( 2 0) and F/ 2 ) for "questionable" links were dropped 



altogether, so one should view this result as a sort of 
"lower bound" for -F(i,o)- Nevertheless, our results are 
a clear indication that the interdimer interactions along 
the (1,0) link in the (ac) plane are asymmetric with re- 
spect to the spin index interchange 1 2. It is worth 
noting that the signs of and i'V-i) agree with the 

pattern proposed in Ref. 

In PHCC, the crystal symmetry makes all J 11 and 
J 22 equal, so the only asymmetry giving rise to nonzero 
amplitudes F$ can be the difference between Jg 2 and 
Jg 1 . The existence of a large Fg for the strongest (1,0) 
link implies that the magnitudes of the individual cou- 
plings Jf^ty are considerably larger than the effective 
"dimer" exchange constant J = Jf-fo) — 0.44 meV. For 
instance, if one assumes for simplicity that ^ = 



1.4J; another, 



then J^ Q) = J 22 Q) ~ 1.7 J and J 21 ^ 
although less likely possibility is that the "diagonal" cou- 
plings have different signs, ^(\ 2 ) < and Jfl \ > 0. 

Magnon decay in S = 1 Haldane chains. - Let us now 
turn to the S — 1 spin chain. Indeed, this model can 
also be viewed as a chain of S = \ dimers with an- 
tiferromagnetic intradimer exchange 4J and an addi- 
tional local constraint symmetrizing the edge spins of 
every two neighboring dimers, and thus coupling them 
into an effective spin-1, see Fig. 01 This constraint 
can be taken into account 'on average' at a mean-field 
level by adding the Lagrange multiplier term of the form 
"A* E n ( S i,»+i ' S 2,n -j), playing the role of a ferromag- 
netic "Hund's rule" interdimer interaction which has to 
be determined self-consistently. It is then obvious that 
the interaction between the dimers is asymmetric and the 
effective Hamiltonian will contain the term of the form 
(jHJ with a single F\ = It is easy to see that this term 
is consistent with all symmetries of the initial Hamil- 
tonian, e.g., it satisfies the reflection symmetry about 
the no-th bond, which in the bond-boson language corre- 



sponds to the transformation t r . 



>— ► — t2n n -r 



as well 



as the symmetry of reflection about the no-th "site" (i.e., 
the interdimer bond to the right of the no-th dimer, see 
Fig- Et) : which corresponds to t n a t— > — t2n +i-n,a- 

Then, a simple golden rule calculation yields the 
magnon linewidth in the Haldane chain as a function of 
the its momentum p: 



r(p) = |M feo , p _ feD | 2 /|4 Q - £ ;_ fe0 |, 



0) 



where ko is the solution of the equation eu +e p ^k = e p> 
and £k — 2(4 J 2 + nJcosfc) 1 / 2 is the single-magnon dis- 
persion. Using the known value of the Haldane gap 
A ~ 0.4J, one can fix the renormalized value of fi to 
be slightly below 4 J. According to JSJ, the matrix ele- 
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a) 



b) 



4J>0 



-u<0 



FIG. 3: a) representation of spin-1 Haldane chain as a spin-i 
chain with symmetrized spin- 1 ; degrees of freedom (denoted 
by the ovals); b) at the mean-field level the symmetrization is 
equivalent to a ferromagnetic "Hund's rule" interaction — fj,. 



ment has the structure Mk, p ~k — fp{k) — f p (j>— k), so 
the linewidth becomes zero at the wave vector threshold 
p = p c , denned by e{p c ) — 2e(p c /2) (here p c € [n, 2ir}); it 
is exactly the point of crossing with the the two-particle 
continuum. Above the continuum threshold, in the vicin- 
ity of p c the linewidth grows as a square root: 



2|/'(p c /2)| 2 e '(p c ) - E '(p c /2) 



|e"(Pc/2)| 



e"(Pc/2) 



1/2 



Vp^Pc- (io) 



It is notable that two-magnon decay processes, strictly 
forbidden in the Haldane chain within the effective con- 
tinuum field theory, nevertheless survive in the full lattice 
description. 

Using the above simple estimate for fi, one can see 
that the magnon linewidth in the Haldane chain grows 
fast to rather large values, e.g., at p = one obtains 
r ~ 3.4J, to be compared to e(p — 0) ~ 5.6J. Such 
a large value of the linewidth arises because the effec- 
tive "interdimer" interaction fi is not small compared to 
the "intradimer" coupling 4 J. An experimental observa- 
tion in that case would show a very fast disappearance of 
quasiparticle excitation above the continuum threshold, 
similar to the picture observed recentl}^ in a spin ladder 



compound (CHa^CHNHaCuCLj. However, one should 
bear in mind that the behavior of linewidth in the vicinity 
of the threshold depends on the details of the wave vector 
dependence of the matrix element, which cannot be ac- 
counted for in the simple qualitative description sketched 
above: for instance, if f'(p c /2) is anomalously small for 
some reason, then one has to take into account higher 
derivatives of f(p) and (|10f) gets replaced by a sum of 
competing terms, T(p) oc X^n=i A n {p — p c ) n ■ Recently, 
a finite magnon linewidth slowly emerging at p c ~ 0.67T 
was observed^ in the S = 1 chain material CSN1CI3 . It 
would be interesting if future experiments could provide 
more information on two-magnon decay in Haldane spin 
chain. 

Summary.- To conclude, we have shown that ele- 
mentary excitations in a gapped dimerized spin sys- 
tem can become unstable above the two-particle contin- 
uum threshold, and that the dependence of the magnon 
linewidth on the wave vector can be used as an additional 
input (complementing the magnon dispersion) to extract 
the information on the exchange interactions pattern. On 
the basis of fitting the recent resultsi on the magnon 
linewidth in the quasi-two-dimensional spin dimer mate- 
rial PHCC, we conclude that interdimer couplings along 
the a axis in this material are strongly asymmetric, con- 
trary to the previously assumed 9 pattern. We argue that 
three-magnon interaction, which makes possible the in- 
stability of a single-magnon excitation against the decay 
into two particles, is generically present in other gapped 
spin systems such as spin-1 Haldane chains. 
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